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ABSTRACT 


Reeent  examinations  of  the  eharaeteristies  of  eoda-derived  S  and  Lg  speetra  for  yield  estimation  have  shown  that  the 
speetral  peak  of  Nevada  Test  Site  (NTS)  explosion  speetra  is  depth-of-burial  dependent,  and  that  this  peak  is  shifted 
to  higher  frequeneies  for  Lop  Nor  explosions  at  the  same  depths.  To  eonfidently  use  eoda-based  yield  formulas,  we 
need  to  understand  and  prediet  eoda  speetral  shape  variations  with  depth,  souree  media,  veloeity  stmeture, 
topography,  and  geologieal  heterogeneity. 

We  are  in  the  final  year  of  a  synthetie  seismogram  study  to  evaluate  the  eharaeteristies  of  Lg  eoda  generated  from 
explosions  and  earthquakes.  To  meet  our  objeetives,  we  have  developed  a  deterministie  model  for  the  western 
United  States  and  the  NTS.  In  order  to  generate  eoda  in  our  2D  and  3D  pseudospeetral  simulations,  heterogeneity 
was  added  in  the  emst  and/or  upper  mantle.  Three  parameters — eorrelation  length,  Hurst  number  (H),  and  fraetional 
veloeity  perturbation  of  the  heterogeneities — are  used  to  eonstruet  different  realizations  of  random  media.  We  have 
estimated  the  heterogeneity  parameters  for  the  NTS  using  available  seismie  and  geologie  data.  Lateral  eorrelation, 
varianee  and  eoherenee  measures  between  seismie  traees  are  estimated  from  elusters  of  nuelear  explosions  and  well- 
log  data.  The  horizontal  von  Karman  eorrelation  length  estimated  for  nuelear  explosion  Lg  phases  is  between 
0.4-1 . 1  km,  with  Hurst  numbers  between  0.3  and  0.6.  A  dominant  heterogeneity  dimension  of  0.5-1 .8  km  results  in 
best  Lg  eoda  eoherenee  at  frequeneies  between  1  and  4  Hz. 

Examples  of  the  Lg  and  Lg  eoda  syntheties  produeed  by  these  models  resemble  NTS  earthquake  and  explosion 
waveforms;  however,  there  is  a  eonsiderable  non-uniqueness  problem  that  must  be  addressed.  We  investigated 
methods  for  quantifying  our  syntheties  in  order  to  assure  that  simulated  waveforms  have  similar  eharaeteristies  as 
observed  eoda  and  direet  Lg.  One  sueh  eharaeteristie  is  that  the  Q  estimates  for  direet  Lg  and  its  eoda  are  typieally 
observed  to  be  near  equivalent  in  many  teetonie  settings.  For  example,  near  the  NTS,  Qig  (direet)  and  Qig  (eoda) 
typieally  range  between  200-300.  Thus,  we  have  eompared  the  Q  values  of  synthesized  Lg  eoda  and  direet  Lg  from 
models  with  different  attenuation  struetures.  These  struetures  range  from  eonstant  shear  wave  Q  (e.g.,  Q^200  or 
500)  throughout  the  erust  to  published,  depth-dependent  models  for  the  Basin  and  Range.  We  examined  the  effeet  of 
the  loeation  of  the  stoehastie  variations  (e.g.,  within  the  entire  erust,  above  the  middle  emst,  or  only  in  the  upper 
emst)  on  the  Q^g  (direet)  and  Q^g  (eoda)  estimates  obtained  from  the  syntheties.  The  eomparison  of  the  variations  in 
Q  estimates  for  different  souree  types  (e.g.,  earthquakes  and  explosions)  at  different  depths  suggests  that  our 
stoehastie  models  produee  realistie  1  Hz  estimates  of  Qig  (direet)  and  Q^g  (eoda)  for  published  Basin  and  Range 
attenuation  models.  Estimates  for  the  frequeney  dependenee  of  the  synthetie  eoda  are  relatively  high  eompared  to 
observed  data. 

We  are  eurrently  using  these  ealibrated  stoehastie  models  to  generate  regional-distanee  syntheties  for  monopole 
explosions  at  depths  ranging  from  0.1  to  1  km  for  all  test  site  models.  We  have  superposed  seeondary  souree  effeets, 
sueh  as  a  eompensated  linear-veetor  dipole  (CLVD)  souree,  on  the  monopole  syntheties.  We  are  now  deriving  Sn 
and  Lg  eoda  speetra  from  the  syntheties,  estimating  moments  and  yields  from  these  speetra,  and  eomparing  them  to 
observed  data  from  NTS.  If  sueeessflil,  this  synthetie  method  may  be  used  to  estimate  the  S  and  Lg  eoda  properties 
for  yield  estimation  of  explosions  at  either  historieal  test  sites  or  for  broad,  unealibrated  regions  where  we  will  likely 
have  little  information  on  veloeity  stmeture. 
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OBJECTIVES 


We  are  attempting  to  determine  why  the  current  mechanisms  for  Sn  and  Lg  generation  from  explosions  often 
produce  stable  coda  magnitudes  that  are  transportable  between  test  sites.  We  aim  to  understand  and  predict  Lg  and 
Sn  coda  spectral  shapes  with  variations  in  source  depth,  material  properties,  velocity  structure,  and  geological 
heterogeneity. 

Our  investigation  objectives  include  the  following: 

1)  Compilation  and  parameterization  of  seismic  velocity  models  consisting  of  deterministic  material  models  with 
stochastic  perturbations; 

2)  Estimation  of  the  effect  of  deterministic  and  stochastic  model  parameters  on  Lg  propagation  using  synthetic 
waveforms  generated  by  the  Generalized  Fourier  Method  (GFM)  (Orrey,  1995); 

3)  Calibration  of  stochastic  variations  at  each  test  site  using  nearby  earthquake  data;  and 

4)  Calculation  of  2D/3D  synthetics  for  composite  explosion  source  models  of  varying  depths  and  moments  for  each 
test  site. 

During  the  past  year,  we  have  worked  on  all  four  of  our  above  objectives  and  are  near  completion  of  the  project.  If 
these  methods  prove  to  be  successful  for  explaining  the  characteristics  of  explosion-generated  direct  Lg  and  Lg  coda 
at  the  NTS,  we  plan  to  test  the  transportability  of  these  modeling  techniques  at  the  Shagan,  Degelen,  Novaya 
Zemlya,  and  Lop  Nor  test  sites. 

RESEARCH  ACCOMPT  JSHED 


Modeling  Direct  Lg  and  Lg  Coda  Q  in  the  Basin  and  Range 

The  objective  of  our  research  is  to  determine  if  numerical  modeling  of  Z.g  and  its  coda  can  replicate  the  depth- 
dependent  spectral  peaking  observed  at  the  NTS  (Murphy  et  ah,  in  review).  In  order  to  complete  this  task,  we  have 
developed  a  velocity  and  attenuation  model  for  the  NTS  and  are  currently  calibrating  the  stochastic  parameters  of 
the  model. 

Velocity  Model.  The  background  deterministic  velocity  structure  for  NTS  (Figure  1)  is  based  on  a  regional  model 
for  the  Basin  and  Range  similar  to  the  model  developed  by  Benz  et  al.  (1991).  The  velocities  in  the  upper  crust  are 
based  on  borehole  data,  geologic  and  gravity  data,  refraction  studies  and  seismic  experiments  (McLaughlin  et  al. 
1983;  Stump  and  Johnson,  1984;  Ferguson  et  al.,  1994;  Stevens  et  al.,  1991).  We  have  used  various  techniques  to 
estimate  stochastic  parameters  for  the  Basin  and  Range.  This  includes  using  previous  well-located  nuclear 
explosions  to  estimate  correlation  lengths  for  scattering  of  the  Lg  phase  (Tibuleac  et  al.,  2005,  2006).  These  results 
suggest  that  the  horizontal  correlation  (Ay)  lengths  for  a  von  Karman  stochastic  model  range  between  0.5  -  1  km 
with  smaller  vertical  correlation  (Az)  lengths.  The  H  numbers  are  related  to  the  fractal  dimension  of  the  medium  and 
range  between  0.3  and  0.6.  Figure  1  shows  the  results  of  converting  the  deterministic  velocity  model  for  the  NTS 
region  into  a  stochastic  model  for  generating  direct  Lg  and  its  coda. 

Attenuation  Models.  We  have  considered  different  intrinsic  background  attenuation  models  in  our  numerical 
modeling  exercises,  ranging  from  constant  shear  wave  Q  models  to  published  Q  models  for  the  western  United 
States.  Bonner  et  al.  (2007)  measured  Qig  (direct)  and  Qig  (coda)  from  synthetics  propagated  through  a  constant 
shear  wave  Q  model  (e.g.,  Qy?  =  200  or  500)  for  the  entire  crust.  Their  modeling  suggested  that  Q^g  (direct)  typically 
ranged  between  the  values  of  Qy?  and  Qa  for  models  without  stochastic  variations.  Stochastic  variations  lowered  the 
apparent  Q^g  (direct);  however,  it  was  troubling  that  the  measured  Qig  (direct)  and  Qxg  (coda)  were  never 
approximately  the  same.  A  comparison  of  regionalized  Q  maps  often  shows  the  near-equivalence  of  these  two 
estimates  of  attenuation. 
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Fignre  1.  Example  of  a  2D  stochastic  model  for  the  NTS  nsed  to  model  Qig  (direct)  and  Qig  (coda)  in  this 
stndy.  Also  Usted  are  the  variations  in  Ay,  Az,  Hnrst  nnmher,  and  pertnrhation  magnitnde  and 
location  nsed  to  make  different  models  for  testing. 


In  this  paper,  we  examine  synthetics  propagated  through  two  published  Q  models  for  the  Basin  and  Range.  The 
Patton  and  Taylor  (PT)  (1984)  model  is  characterized  by  low  (85-172)  throughout  the  entire  crust.  Similarly,  the 
Mitchell  and  Xie  (MX)  (1994)  model  has  low  in  the  upper  crust  (50-80);  however,  they  suggest  high  in  the 
middle  to  lower  cmsts  (400-1000).  We  incorporated  these  two  Q  models  into  our  deterministic  and  stochastic 
models  (Figure  1)  for  the  NTS  and  computed  synthetics  for  earthquakes  using  the  GFM  (Orrey,  1995). 

Methodology.  Record  sections  of  GFM  synthetic  waveforms  generated  by  four  different  perturbations  of  our  NTS 
model  are  shown  in  Figure  2,  including  deterministic  (e.g.,  without  stochastic  heterogeneity)  and  stochastic  (Ax  =  1, 
Az  =  0.5)  versions  for  both  the  PT  and  MX  attenuation  models.  The  appearance  of  the  synthetics  with  stochastic 
variations  compares  favorably  with  observed  waveforms  from  earthquakes  and  explosions  in  the  Basin  and  Range 
(see  Tibuleac  et  al.,  2005,  2006  for  additional  examples). 

Direct  Lg  Q  Estimation.  We  windowed  the  direct  Lg  arrivals  between  3.6  and  3.0  km/sec  and  estimated  the  spectral 
amplitudes  A  at  each  station  at  distance  R.  We  used  a  cosine  taper  and  smoothed  the  spectra  in  order  to  reduce  the 
variance  between  the  stations.  For  each  frequency  / we  used  the  following  equation  (Mitra  et  al.,  2006): 

logio  A  +  Q-5  logio  7?  ^  logio  5  -  ^  ^  (2) 

to  solve  for  Q(f).  For  each  frequency,  we  plotted  logioT  +  0.5  log  lo  R  versus  R  and  performed  a  linear  regression  to 
determine  Q.  We  used  a  velocity  (v)  of  3.5  km/sec  for  the  direct  Lg  arrival.  Plots  of  the  Q  as  a  function  of  frequency 
are  then  used  to  estimate  the  Q  at  1  Hz  (Qo)  and  the  frequency  dependence  (rj)  (see  Figure  3).  Because  of  the  grid 
size  and  velocities  in  our  models,  our  simulations  are  only  valid  at  frequencies  less  than  2-3  Hz. 
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Figure  2.  GFM  Synthetic  waveforms  filtered  near  1  Hz  for  deterministic  (Top  Row)  and  stochastic  (Bottom 
Row)  models  for  the  NTS.  Synthetics  in  the  left  column  were  propagated  through  the  Mitchell  and 
Xie  (1994)  Q  model  while  the  Patton  and  Taylor  (1984)  Q  model  was  used  to  attenuate  the  synthetics 
in  the  right  column. 


Ay=1  Az=0.5  20%  with  MX  Q  Model  Ay=1  Az=0.5  20%  with  PT  Q  Model 


Figure  3.  Estimating  Qig  (direct)  from  synthetics.  (Left)  Regression  of  log  Q  versus  frequency  for  Lg 

propagated  through  a  stochastic  NTS  model  with  the  Mitchell  and  Xie  (1994)  shear-wave  Q  model. 
(Right)  Same  as  left  except  for  use  of  Patton  and  Taylor  (1984)  Q  model.  There  appears  to  he  two 
different  trends  (above  and  helow  1  Hz);  however,  for  this  study,  we  can  fit  most  of  the  data  plus 
error  with  a  single  line. 
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Given  the  non-uniqueness  associated  with  synthetic  modeling  efforts  similar  to  this,  we  decided  to  determine  if  our 
GFM  synthetics  have  similar  characteristics  as  observed  coda  and  direct  Lg.  In  the  remaining  sections,  we  compute 
the  apparent  (direct)  and  Qig  (coda)  values  for  the  synthetics  using  different  Q  models  for  comparison  to 
observed  estimates  in  the  Basin  and  Range. 

Lg  Coda  Q  Estimation.  To  estimate  Qig  (coda) ,  we  used  the  stacked  spectral  ratio  (SSR)  method  of  Xie  and  Nuttli 
(1988).  In  this  method,  amplitude  spectra  for  windows  at  different  times  within  the  synthetic  coda  are  used  to  form 
spectral  ratios.  The  ratios  are  corrected  for  velocity  and  geometrical  spreading,  stacked,  and  regressed  versus 
frequency  to  determine  the  coda  Q  at  1  Hz  (Q^)  and  the  frequency  dependence  {rf}.  Examples  of  the  SSR  analysis  on 
synthetic  data  are  shown  in  Figure  4. 

Results.  Table  1  provides  measured  (direct)  and  Qig  (coda)  in  the  Basin  and  Range  from  previous  studies.  These 
results  suggest  that  both  the  direct  and  coda  Q  for  the  Lg  phase  recorded  in  this  region  should  range  between  200- 
300  with  T]  between  0.40  and  0.80.  Table  2  and  Figure  5  show  the  results  of  measuring  the  direct  Lg  Q  and  Lg  coda 
Q  for  synthetics  propagated  through  our  different  deterministic  and  stochastic  models  for  the  Basin  and  Range. 
Without  the  presence  of  stochastic  variations,  GFM  synthetics  in  an  XM  Q  model  do  not  produce  Qxg  (direct)  values 
that  match  the  measured  estimates  for  the  Basin  and  Range.  However,  with  20%  stochastic  variations  added  to  the 
model,  the  XM  Q  model  produces  synthetics  with  Q^g  (direct)  and  Q^g  (coda)  that  fall  within  the  range  of  observed 
data.  The  synthetic  Q^gS  (direct)  for  all  PT  models  were  within  the  measured  range;  however,  adding  20%  variations 
resulted  in  Q^g  (coda)  estimates  that  were  slightly  below  observed.  Unfortunately,  our  synthetics  do  not  reproduce 
the  observed  Q^g  (coda)  frequency  dependence  rj  no  matter  which  Q  model  is  used.  In  the  coming  months,  we  hope 
to  determine  the  reason  for  this  discrepancy,  which  could  include  the  effects  of  using  2D  versus  3D  models  as  well 
as  possible  need  for  frequency-dependent  intrinsic  attenuation  in  the  model. 


Ay=1  Az=0.5  20%  with  MX  Q  Model  Ay=1  Az=0.5  20%  with  PT  Q  Model 


Frequency  (Hz)  Frequency  (Hz) 

Figure  4.  SSR  method  (Xie  and  Nuttli,  1988)  applied  to  synthetic  waveform  data  for  Lg  coda  Q  estimation  for 
a  stochastic  model  with  Mitchell  and  Xie  (left)  and  Patton  and  Taylor  attenuation  (right). 
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Table  1.  Ohserved  (direct)  and  (coda)  for  the  Basin  and  Range 


Q  Study 

Qig  (direct) 

Qig  (coda) 

Lg  Qo  (direct) 

n 

Singh  and  Herrmann  (1983) 

200-300 

0.40-0.60 

Chavez  and  Priestley  (1986) 

214±  15 

0.54 

Xie  and  Mitehell  (1990) 

267  ±56 

275  ±  26 

Baqer  and  Mitehell  (1998) 

250-300 

0.40-0.60 

Benz  et  al.  (1997) 

235  ±  11 

0.56  ±0.04 

Aleqabi  and  Wysession  (2006) 

234-312 

0.40-0.80 

Table  2,  Synthetic  Qls  (direct)  and  (coda)  for  Basin  and  Range  models 


Qig  (direct)  (coda)  Stochastic  Parameters 


MODEL 

Lg  Qo 

Error 

0 

Error 

CODA 

Qo 

Error 

0 

Error 

Q  Model 

Ay 

Az 

Pereent 

1 

225 

22 

1 

0.05 

214 

13 

1.1 

0.1 

PT 

1 

0.5 

10% 

2 

207 

23 

0.81 

0.09 

195 

12 

0.92 

0.09 

PT 

1 

0.5 

20% 

3 

253 

26 

0.9 

0.05 

208 

16 

1.16 

0.1 

PT 

1 

0.25 

10% 

4 

207 

21 

0.96 

0.03 

197 

9 

0.88 

0.08 

PT 

1 

0.25 

20% 

5 

267 

22 

1.3 

0.21 

... 

- 

— 

PT 

... 

— 

... 

6 

385 

51 

0.69 

0.04 

263 

25 

1.11 

0.11 

XM 

1 

0.5 

10% 

7 

314 

41 

0.57 

0.07 

237 

20 

0.84 

0.11 

XM 

1 

0.5 

20% 

8 

370 

59 

0.98 

0.07 

251 

20 

1.17 

0.1 

XM 

1 

0.25 

10% 

9 

310 

31 

1.1 

0.06 

233 

13 

0.91 

0.1 

XM 

1 

0.25 

20% 

10 

632 

127 

1 

0.05 

... 

- 

— 

XM 

... 

— 

... 

Q  at  1  Hz  from  Synthetics 
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Fignre  5.  Qo  (left)  and  t)  (right)  determined  for  Lg  synthetics  from  Basin  and  Range  stochastic  models 
compared  to  observed  estimates  (shaded  regions). 

Modeling  NTS  Explosion-Generated  Coda.  The  earthquake  simulations  using  the  Patton  and  Taylor  (1994) 
Q  model  for  the  Basin  and  Range  together  with  various  stoehastie  perturbations  provide  Qig  (direet  and  eoda) 
estimates  that  best  mateh  observed  data  (Table  2).  We  generated  explosion  syntheties  for  PT  models  1-4  and 
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compared  them  to  NTS  nuclear  explosions  recorded  at  station  MNV  (distance  ~230  km).  Figure  6  shows  a 
comparison  of  the  MNV  recordings  of  the  NTS  event  Brie  compared  with  synthetics  for  Model  1.  Model  1  consists 
of  the  PT  Q  model  and  stochastic  variations  up  to  10%  with  correlation  distances  of  1.0  to  0.5  km.  The  synthetics 
and  observed  data  were  normalized  to  their  maximum  Pg  amplitudes.  Brie  was  detonated  a  depth  of  less  than  250 
meters  (Murphy  et  ah,  in  review);  therefore,  we  placed  our  monopole  explosion  source  at  a  depth  of  250  meters  in 
order  to  compute  the  synthetics.  The  synthetics  are  valid  at  frequencies  less  than  2-3  Flz,  and  both  the  observed  and 
synthetic  waveforms  in  Figure  6  have  been  bandpass  fdtered  between  0. 1  and  2  Flz. 

A  quick  glance  of  the  waveforms  in  Figure  6a  might  suggest  good  agreement  between  the  observed  and  synthetic 
arrivals.  The  Pn  and  Pg  arrivals  show  very  good  correlation  in  both  relative  amplitude  and  arrival  times.  We  note 
that  the  Lg  arrival  (green  box)  is  not  as  pronounced  in  the  synthetic  waveform  as  it  is  in  the  observed  data,  and  that 
the  synthetic  coda  (blue  box)  has  higher  frequency  content  and  more  rapid  amplitude  decay  than  the  observed  data. 
This  is  also  shown  by  the  coda  spectra  in  Figure  6b  which  shows  a  large  amplitude  peak  near  0.9  Flz  in  the  observed 
data  that  is  missing  from  the  synthetic  coda. 


a)  1 


Figure  6.  Explosion  synthetics  and  data.  a).  Comparison  of  observed  waveforms  for  the  NTS  explosion  Brie 
recorded  at  NTS  compared  to  GEM  synthetics  for  the  explosion  using  Model  1  from  Table  2. 

The  time  windows  for  Lg  and  Lg  coda  (green  and  blue  boxes,  respectively)  were  used  to  compute 
spectra  (b). 
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The  correlation  distances  for  the  Model  1  perturbations  were  1  and  0.5  km,  and  the  results  shown  in  Figure  6  suggest 
that  coda  energy  at  frequencies  less  than  1  Flz  was  being  scattered  too  rapidly  when  compared  to  observed  data.  By 
decreasing  the  correlation  distance  from  0.5  to  0.25  km,  we  should  decrease  the  rapid  scattering  below  1  Flz.  This  is 
confirmed  in  the  seismograms  and  spectra  for  simulations  using  Model  4  (Figure  7).  The  simulations  for  Model  4 
have  increased  energy  near  1  Flz  in  the  first  25  seconds  of  the  coda  analysis  window  compared  to  the  synthetics 
shown  in  Figure  6.  We  note  peaks  in  the  synthetic  and  observed  spectra  between  0.9  and  1  Flz  (Figure  7b).  These 
synthetics  provide  an  excellent  comparison  to  the  observed  data  for  the  P,  P  coda,  direct  Lg,  and  early  coda  arrivals. 
Our  synthetics  do  not  do  an  adequate  job  of  modeling  the  observed  coda  at  periods  less  than  0.7  Flz  at  onset  times 
greater  than  120  seconds  (see  Figure  8),  which  may  be  the  result  of  reaching  the  limit  of  our  2D  coda  modeling 
methods.  The  later  arriving  coda  may  be  multi-pathing  from  structures  away  from  the  direct  path,  which  can  not  be 
effectively  modeled  using  only  2D  profiles.  In  the  remaining  months  of  the  project,  we  will  also  add  a  CLVD  source 
to  the  monopole  explosion  to  see  if  the  addition  of  more  short-period  surface-wave  energy  could  improve  the  match 
between  observed  and  synthetic  coda. 


0  50  ,  100  150 

Time  (sec) 


Figure  7.  Explosion  synthetics  and  data,  a)  Comparison  of  observed  waveforms  for  the  NTS  explosion  Brie 
recorded  at  NTS  compared  to  GEM  synthetics  for  the  explosion  using  Model  4  from  Table  2. 

The  time  windows  for  Lg  and  Lg  coda  (green  and  blue  boxes,  respectively)  were  used  to  compute 
spectra  (b). 
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Synthetic  Brie 


Seconds  Seconds 

Fignre  8.  A  comparison  of  the  coda  envelopes  for  the  Model  4  synthetics  (left)  and  observed  explosion 
seismograms  (right)  for  five  narrow  hands  between  0.3  to  2  Hz. 

CONCLUSIONS  AND  RECOMMENDATIONS 

Our  objective  is  to  try  to  develop  synthetic  seismograms  for  5-wave  coda  from  explosions.  The  ability  to  accurately 
model  direct  Lg  and  Lg  coda  using  accurate  earth  models  (both  velocity  and  attenuation)  will  be  important  for 
regions  in  which  explosion  calibration  data  do  not  exist.  Our  results  for  NTS  suggest  we  can  numerically  model 
direct  Lg  and  Lg  coda  that  have  similar  attenuation  and  propagation  characteristics  as  observed  data.  We  can  do  so 
using  literature-based  velocity  and  attenuation  models  and  standard  estimates  for  correlation  lengths.  We  also  have 
modeled  explosion-generated  coda  that  have  attributes  such  as  peak  frequency  and  relative  amplitudes  that 
sometimes  closely  match  observed  nuclear  explosion  data,  especially  for  the  early  onset  arrivals  in  the  coda 
windows. 

Currently,  our  models  do  not  produce  realistic  estimates  of  the  frequency  dependence  for  Lg  coda  in  the  Basin  and 
Range.  Aki  (1969)  described  the  coda  as  the  result  of  backscattered  waves  from  numerous  heterogeneities 
distributed  randomly  and  uniformly  in  the  crust  and  upper  mantle  within  an  ellipse.  Based  on  his  definition,  the 
scattering  that  creates  the  coda  is  from  3D  phenomena;  this  suggests  our  initial  modeling  in  2D  may  never  provide 
us  with  a  perfect  match  between  observed  and  synthetic  coda.  In  order  to  achieve  a  model  where  the  synthetic  coda 
matches  the  observed,  we  may  have  to  consider  3D  models;  however,  this  will  ultimately  either  increase  the 
computation  time  greatly  or  decrease  the  resolvable  frequency  bandwidth  to  below  1  Hz.  Example  synthetics  for  Lg 
and  Lg  coda  from  our  first  attempts  at  incorporating  3D  stochastic  models  are  shown  in  Figure  9. 

3D  STRIP  MODEL22  BP;  0.8  to  1.2Hz 


Figure  9.  GEM  synthetics  in  a  3D  stochastic  velocity  model.  The  Lg  group  velocity  window  is  outlined  using 
the  red  lines. 
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